The thawing dark energy dynamics: Can we detect it? 
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We consider different classes of scalar field models including quintessence and tachyon scalar fields 
with a variety of generic potential belonging to the thawing type. We focus on observational quanti- 
ties like Hubble parameter, luminosity distance as well as quantities related to the Baryon Acoustic 
Oscillation measurement. Our study shows that with present state of observations, one can not 
distinguish amongst various models which in turn can not be distinguished from cosmological con- 
stant. Our analysis indicates that there is a thin chance to observe the dark energy metamorphosis 
in near future. 

PACS numbers: 98.80 Cq 



I. INTRODUCTION 

The fact that our universe is currently going through 
an accelerated phase of expansion is one of the most sig- 
nificant discoveries [l[ in physics in recent times that can 
have far reaching implications for fundamental theories 
of physics. Late time cosmic acceleration can be fueled 
cither by assuming the presence of an exotic fluid with 
large negative pressure known as dark energy or by mod- 
ifying gravity itself. The simplest candidate of dark en- 
ergy is provided by cosmological constant with equation 
of state parameter w = — 1. However, the model based 
upon cosmological constant is plagued with the fine tun- 
ing and cosmic coincidence problems (See Ref Q for a 
nice review). 

Scalar field models with generic features can alleviate 
the fine tuning and coincidence problems and provide 
an interesting alternative to cosmological constant Q. 
The simplest generalization of cosmological constant is 
provided by a scalar field with linear potential Its 
evolution begins from the locking regime (due to large 
Hubble damping) where it mimics the cosmological con- 
stant like behavior. At late times, the field starts rolling 
and since the potential has no minimum, the model leads 
to a collapsing universe with a finite history. 

The more complicated scalar field models can broadly 
be classified into two categories- fast roll and slow roll 
models dubbed freezing and thawing models In case 
of the fast roll models, the potential is steep allowing the 
scalar field to mimic the background being subdominant 
for most of the evolution history. Only at late times, the 
field becomes dominant and drives the acceleration of the 
universe. Such solutions are referred to as trackers. 

Slow-roll models are those for which the field kinetic 
energy is much smaller than its potential energy. It usu- 
ally has a sufficiently flat potential similar to an inflaton. 
At early times, the field is nearly frozen at w = — 1 due 
to the large Hubble damping. Its energy density is nearly 
constant and and its contribution to the total energy den- 
sity of the universe is also nearly negligible. But as ra- 
diation/matter rapidly dilutes due to the expansion of 
the universe and the background energy density becomes 
comparable to field energy density, the field breaks away 



from its frozen state evolving slowly to the region with 
larger values of equation of state parameter. In this case, 
however, one needs to have some degree of fine tuning 
of the initial conditions in order to achieve a viable late 
time evolution. 



Recent observations suggest that the equation of state 
parameter for dark energy does not significantly deviate 
from w = — 1 around the present epoch @. This type 
of equation of state can be easily obtained in dynam- 
ical models represented by thawing scalar fields. This 
fact was exploited in Ref |7( which examined quintessence 
models with nearly flat potentials satisfying the slow-roll 
conditions. It was shown that under the slow-roll con- 
ditions, a scalar field with a variety of potentials V{(j>) 
evolve in a similar fashion and one can derive a generic 
expression for equation of state for all such scalar fields. 
This result was later extended to the case of phantom 
Q and tachyon scalar fields Q. It was demonstrated 
that under slow-roll conditions, all of them have identi- 
cal equation of state and hence can not be distinguished, 
at least, at the level of background cosmology. The cru- 
cial assumption, for arriving at this important conclusion 
was the fulfillment of the slow-roll condition for the field 
potentials. 

In this paper, we relax the assumption of slow roll but 
assume that the scalar field is of thawing type i.e it is 
initially frozen at w = —1 due to large Hubble damp- 
ing. With non-negligible matter contribution, this does 
not necessarily mean the small value for V^/V which is 
the usual slow-roll parameter for inflaton. We, rather, 
assume that the slow-roll condition is highly broken such 
that V^/V ~ 1. In this case, we need to fine tune the 
initial conditions to match the observational value of the 
present day dark energy density which is a characteris- 
tic feature of any thawing model. With this choices, we 
study the evolution of a variety of scalar field models 
having both canonical and non-canonical kinetic terms. 
We particularly focus on the observational quantities like 
Hubble parameter, luminosity distance as well as quan- 
tities related to the Baryon Acoustic Oscillation (BAO) 
measurement. 



2 



II. THAWING SCALAR FIELD 



where 



In what follows, we shall assume that the dark energy 
is described by a minimally-coupled scalar field, 0, with 
equation of motion 

4> + 3H0 + V^ = Q (1) 

where the Hubble parameter H is given by 



(2) 



Here p is the total energy density in the universe. We 
model a flat universe containing only matter and a scalar 
field, so that 0^, + Ojvf = 1. 

Equation (Q]) indicates that the field rolls downhill in 
the potential V((f>), but its motion is damped by a term 
proportional to H. The equation of state parameter w 
is given by w = p$/ Prf, where the pressure and density of 
the scalar field have the form 



(3) 
(4) 



Observations suggest a value of w near —1 around the 
present epoch. We adopt a similar technique as followed 
in references 0, H| and define the variables x, y, and A as 



y = ^v{4>)/m\ 

1 dV 



(5) 

(0) 
(7) 



where prime as usual denote the derivative with respect 
to In a: e.g., <f>' = a(d(f)/da) 

Then contribution of the kinetic energy and potential 
energy of the scalar field to the fractional density param- 
eter D,^ are represented by x 2 and y 2 such that, 



VL^ = x 2 + y 2 , 
while the equation of state is given by, 

2x 2 



7 = 1 + w 



x 2 + y 2 



(8) 



(9) 



In terms of the variables x, y, and A, evolution equa- 
tions ([T]) and (J2J> take the autonomous form 



-3x + \J'^y 2 + ^x[l + x 2 -y 2 }, 



y = -Ay -xy + -y[i + : 
A' = -V6A 2 (r-I)x, 



(10) 

(11) 
(12) 



dej) 2 ' V# 



dV 



(13) 



We now rewrite these equations, changing the dependent 
variables from x and y to the observable quantities fi^ 
and 7 given by equations ((8]) and (J9j) . To make this trans- 
formation, we assume that x' > 0; our results generalize 
trivially to the opposite case. In terms of Q<p and 7 the 
above set of equation become 



7' = -37(2-7) + A(2- 7 )v/37^, 
% = 3(1- 7 )fy(l -fy), 

-V3A 2 (r- i)Vt^V 



A' 



(14) 
(15) 

(16) 



This is an autonomous system of equations involving 
the observable parameter 7 and fi^. Given the initial 
conditions for 7, and A, one can solve this system 
of equation numerically for different potentials. As we 
mention earlier, we are interested in thawing models i.c 
models for which the equation of state is initially frozen 
at w = —1. Hence 7 = initially for our purpose. We 
also do not assume slow-roll conditions for the scalar field 
potentials rather we consider situations for which it is 
broken strongly i.e Xinitial ~ 1- We should mention that 
for models where slow- roll condition is satisfied i.e. A << 
1, it has been already shown that all such models have 
an identical equation of state as a function of scale factor 
0, Hi ■ In general the contribution of the scalar field to 
the total energy density of the universe is insignificant at 
early times, nevertheless one has to fine tune the initial 
value of f2<£ in order to have its correct contribution at 
present. This is the fine tuning one needs to have in a 
thawing models. With these initial conditions we evolve 
the above system of equations from redshift z = 1000 
(or a = 10~ 3 ) till the present day z = (a = 1). We 
consider various types of potentials e.g V = <j), V = (f> 2 , 
V = and V = (f>~ 2 , characterized by T = 0, 5, 1 and 
I respectively. We have taken two sets of solution such 
that ft^ = 0.7 and il* = 0.75 at the present epoch for all 
chosen values of T. 

We also consider the Pseudo-Nambu Goldstone Boson 
(PNGB) model [lo|. (For a recent discussion, see Ref. 
[llj and references therein). This model is characterized 
by the potential 



V(0) = M 4 [cos(0//) + l], 



(17) 



Alam et al.,[12j have previously considered such type of 
potential to see whether dark energy is decaying or not. 

We have chosen / to be 1 for our purpose without any 
loss generality. 

As mentioned before, tachyon field is of interest in cos- 
mology. There have been several investigations using this 
field as a dark energy candidate (3]. In what follows, 
we shall repeat the above presented analysis for tachyon 
field. 
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FIG. 1: Plot of equation of state w vs. redshift for different 
scalar field and tachyon models. Solid curves represent differ- 
ent Tachyonic models with V(<j)) = <j), (j) 2 , e 1 *, cf>~ 2 respectively 
from top to bottom, Dashed curves from top to bottom rep- 
resent different scalar field models with same potentials as in 
tachyon. Dotted curve represents PNGB model. Q m = 0.3. 



III. THAWING TACHYON FIELD 

The tachyon field is specified by the Dirac-Born-Infeld 
(DBI) type of action [13| 



S = / -V(4>)y/1 - d^W^y/^gcPx. 



(18) 



In FRW background, the pressure and energy density of 
the tachyon field <f> are given by 



V{4>) 



(19) 



(20) 



1 



The equation of motion which follows from (|18[) is 

4> + 3H<j>(l-<p 2 ) + ^-(l- l j ) 2 ) = (21) 



where H is the Hubble parameter. The evolution equa- 
tions can be cast in the following autonomous form for 
the convenient use 



x' t = -(1 - x 2 ){3x t - VS\ t y t ) 



, Vt 
Vt= 2 



-Vsx 



txtyt 



X't = -V3X 2 x t y t (T 



(22) 
(23) 

(24) 



FIG. 2: Plot of h(a) = 2<&L V s. redshift for different models 
as well as for ACDM corresponding to f2 m o = 0.3. 



with xt,yt,Xt and T defined as 



Xt = 4>, y t = ^f^ a 4 



(25) 



where prime again denotes the derivative with respect to 
ln(a). We would further use the following definitions for 
the tachyon field as we did in case of thawing quintessence 



fid = 



V 



7=1 



(26) 



where w — — is the equation of state for the tachyon 

Pip 

field. One can now express the autonomous equations 
through them: 



-6 7 (1 



■7)+ 2 v /3 7 fi A t (1-7)2 
ffy =3fid(l-7)(l-^) 

Ai = -V37n;A a (i-7)4(r-|) 



(27) 
(28) 

(29) 



We adopt a similar treatment to solve the above set 
of equation (f2"7} - ([2l?)) as we had done in the earlier case. 
Infact, we even consider similar kind of potentials for 
tachyon fields as well, i.e, V = cf>, V = (j> 2 , V = and 
V = <fi~ 2 , characterized by r = 0, |, 1 and | respectively 
along with the initial conditions for Xunitial to be 1 and 
^initial ~ 0. Here also we take two solutions set for all 
r's, with two different initial conditions of fi^ such that 
at present it contributes 70% and 75% of the total energy 
share. Before discussing our result, we want to point 
that the system of equations (14)-(16) and (27)-(29) for 
scalar field and tachyon field respectively are completely 
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FIG. 3: Plot of A/i(explained in the text) vs. redshift for 
different models. Solid curves represent different Tachyonic 
models with V((f>) = cj>, cj> 2 , e , respectively from top to 
bottom, Dashed curves from top to bottom represent different 
scalar field models with same potentials as in tachyon. Dotted 
curve represents PNGB model. il m o = 0.3. The black dots 
represents the value of the error bar as quoted in [l5|, . We 
assume ho = 73 ± 3kms~ Mpc" 1 for our purpose. 



different. Hence a priori one expects to have different 
evolutions for different potentials as well as for scalar 
and tachyon fields. 

Once wc know the solution for fi^(a) by solving either 
(14)-(16) or (27)-(29), we can easily find the behavior 
of the Hubble parameter which, in terms of Ct^, can be 
expressed as 



h 2 (a) = 



H 2 (a) l-fl 



Hi 



i - n 



(30) 



where Hq and Q^o are the present day values for the 
Hubble parameter and the dark energy density parame- 
ter. This is the most important parameter as all the ob- 
servable quantities involving background cosmology can 
be constructed from this. Moreover there arc indepen- 
dent observational constraint on this parameter itself. 
We should mention that in this approach, one does not 
need to know the equation of state (7(a) = 1 + w(a)) to 
construct the observational quantities although its effect 
comes through the solutions of Eqs.(14)-(16) or (27)-(29). 



IV. RESULTS 

Let us now discuss the results of our investigations. 
In Figl, we plot the behavior of the equation of state 
w(a) for different models. It shows that the equation of 



FIG. 4: Plot of Ah vs. redshift for different models same as 
Fig 3 but with Si m0 = 0.25. 



states of different fields with different potentials behave 
differently as one approaches the present day although 
in the past their behavior are almost identical which is 
not surprising as we have assumed the violation of slow- 
roll condition, i.e Xmitiai ~ 1- With slow- roll condition 
satisfied, i.e, << 1, it was shown earlier that models 
with different potentials have the identical w(a) both for 
scalar and tachyon fields 0-Q ■ 

Next we investigated the behavior of the Hubble pa- 
rameter in different cases. With different behavior for 
w(a) as shown in Figl, one would expect to have different 
behavior for Hubble parameter also. However, they are 
completely indistinguishable as shown in Fig 2. In this 
figure, we have also plotted h(a) for ACDM. It is seen 
that we can not differentiate between individual mod- 
els as well as these models from ACDM. It is interesting 
to note that despite having completely different behavior 
for equation of state, all the models have almost identical 
evolution for the Hubble parameter. This is crucial as all 
the observational quantities are constructed out of h(a) 
at least at the level of background cosmology. In recent 
past, estimates of H(z) were derived by Simon, Verde 
and Jimenez using passively evolving galaxies [15| (also 
see Ref[16]). Keeping this in mind, we next plot the dif- 
ferent Ah = hfieid — hACDM for each of our model in Fig. 
3 and Fig. 4 for two different values of f2 m o (0.3 and 0.25) 
respectively. In the same figures, we have also plotted the 
values of the error bars Ah (where Ah = 



HAH n 
T^- 



We take the values of H and AH from the the data 
(l5l flGj and use the prior H = 73 ± 3kms~ 1 Alpc~ 1 as 
quoted therein. Fig. 3 &4 show that with Cl m0 = 0.3, 
one can not distinguish all the models from ACDM as 
the difference between them is much smaller than the 
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FIG. 5: Plot of A/x(explained in the text) vs. redshift for 
different models. Solid curves represent different Tachyonic 
models with V(cj>) = 4>, 4> 2 , e , <fi~ 2 respectively from top to 
bottom, Dashed curves from top to bottom represent different 
scalar field models with same potentials as in tachyon. Dot- 
ted curve represents PNGB model. Q m o = 0.3. The black 
dots represent the value of the error bars as quoted in the 
Constitution data set 11711 - 



present error bars. With smaller values of the density 
parameter, f2 m o = 0.25, the difference becomes larger 
but the error bars still do not allow to distinguish the 
models from ACDM. The other interesting feature which 
one notices from both these figures, is that the difference 
is maximum in the redshift range from z = 0.5toz = l. 
Hence having more data points for higher redshifts may 
not be useful for the purpose to distinguish different mod- 
els from ACDM; rather the low redshift measurements is 
more vital to do the needful. 



Next we consider the Supernova Type-la observation 
which is one of the direct probes for late time accelera- 
tion. It measures the apparent brightness of the super- 
novae as observed by us which is related to the luminosity 
distance (1l(z) defined as 



(1 + *) 



dz' 
o W)' 



(31) 



FIG. 6: Plot of Afj, vs. redshift for different models same as 
Fig 5 corresponding to Q m o = 0.25. 
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Table 1 

With this, one constructs the distance modulus fi 
which is experimentally measured: 



a = m — M = blog 



d L 
Mpc 



25, 



(32) 



where m and M are the apparent and absolute magni- 
tudes of the supernovae which are logarithmic measure 
of flux and luminosity respectively. In Fig 5 and Fig 6, 
we plot the difference Afi = (1$ — fJ>\c dm for each model 
together with current error bars as quoted in the latest 
Constitution data set (l7j . One can now sec that with 
fl m0 = 0.3 (Fig 5), the difference with ACDM for any 
model is quite small as compared to the value of the er- 
ror bars. But with O m o = 0.25, this difference enhances, 
and models like tachyon and scalar field with linear po- 
tentials as well as scalar field with PNGB potential, have 
significant difference with ACDM which is in the range 
of the values of the error bars. The plots also shows that 
the intermediate redhsift range between 0.4 and 1.0 is 
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most sensitive to compare the models with ACDM and 
future experiments involving Type-la supernova should 
focus more in this rcdshift range in order to investigate 
the nature of dark energy. 
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Table 2 

Another observational probe that has been widely used 
in recent times to constraint dark energy models is re- 
lated to the data from the Baryon Acoustic Oscillations 
measurements 181. In this case, 



one needs to calculate 
the parameter Dy which is related to the angular diam- 
eter distance as follows 



Dv(zbao) 



zbao 
H(zbao) 



dz 
H(zj 



1/3 



(33) 



For BAO measurements we calculate the ratio Dy(z = 
0.35) /LV(z = 0.20). As shown in HU this ratio is a 
relatively model independent quantity and has a value 
1.812 ± 0.060. For our case, we calculate the difference 
of this ratio between any scalar field model and ACDM 
model. In tables 1 & 2, we quote our result for two 
different values for Q m0 : Cl m = 0.3, Q m o = 0.25 . 

As one can see from the results quoted in these two 
tables, with current BAO measurements, it is hard to 
distinguish all the models from ACDM. One has to de- 
crease the error bars at least by fifty percent for this 
purpose. 



V. CONCLUSION 

In this paper, we have studied the general classes of 
thawing models with both quintessence and tachyon type 



scalar fields without assuming slow-roll conditions for the 
potentials of these fields. Our investigations show that 
the overall Hubble parameter has almost identical behav- 
ior for all these models and also matches with its coun- 
terpart corresponding to ACDM despite of the fact that 
the equation of state for different models behaves quite 
differently. Since all the observable quantities related to 
background evolution, are constructed out of H(z), it is 
practically impossible to distinguish these models from 
ACDM using the current data. While analyzing the ob- 
servational constraints, we used supernova and BAO data 
along with the information on Hubble parameter mea- 
surements. Our analysis shows that for smaller values 
of fl m o, it is easier to distinguish amongst the various 
models. We find that tachyon and scalar field models 
with linear potential as well as scalar field model with 
PNGB potential are comparatively easier to distinguish 
from ACDM. It should be mentioned that little attention 
is paid in the literature to scalar models with linear po- 
tentials [20l| : these models deserve further investigations. 

An interesting outcome of our study is related to the 
small redshift range showing distinguished features. We 
find that the behavior of thawing dynamics around the 
rcdshift interval, z = 0.5 —1, is most sensitive to study 
the deviations from ACDM. In our opinion, future obser- 
vations should concentrate more on this particular range 
so as to decrease the error bars significantly. 



VI. ACKNOWLEDGEMENT 



A.A.S acknowledges the financial support provided 
by the University Grants Commission, Govt. Of In- 
dia, through major research project grant (Grant No:33- 
28/2007(SR)). S.S acknowledges the financial support 
provided by University Grants Commission, Govt. Of 
India through the D.S.Kothari Post Doctoral Fellowship. 



[1] R. A. Knop et al., Ap.J. 598, 102 (2003); A. G. Riess, et 
al., Ap.J. 607, 665 (2004). 

[2] E. J. Copeland, M. Sami and S. Tsujikawa, 
Int.J.Mod.Phys.D 15, 1753 (2006); M. Sami, 
larXiv:0904.3445l V. Sahni and A. A. Starobinsky, 
Int.J.Mod.Phys.D 9, 373 (2000); T. Padmanabhan, 
Phys.Rep. 380, 235 (2003); E. V. Linder, astro- 



ph/0704.2064; J . Frieman, M. Turner and D. Huterer, 
arXiv:0803.0982 R. Caldwell and M. Kamionkowski, 
larXiv:0 903.0866 A. Silvestri and M. Trodden, 
larXiv:0904.0024l 
[3] B. Ratra and P. J. E. Peebles, Phys.Rev.D 37, 3406 
(1988); R. R. Caldwell, R. Dave and P. .J .Stein- 
hardt, Phys.Rev.Lett. 80, 1582 (1998); A. R. Liddle 



7 



and R. J. Scherrer, Phys.Rev.D 59, 023509 (1999); 

P. J. Steinhardt. L. .Wang and I. Zlatev, Phys.Rev.D 

59, 123504 (1999) 
[4] R. Kallosh, et al., JCAP 0310, 015 (2003); P. P. Avelino, 

et al., Phys.Rev.D 70, 083506 (2004). 
[5] R. R. Caldwell and E. V. Linder, Phys.Rev.Lett. 95, 

141301 (2007). 

[6] W. M. Vassey, et al., Ap.J. 666, 694 (2007); T. M. Davis, 

et al., Ap.J 666, 716 (2007). 
[7] R. J. Scherrer and A. A. Sen, Phys.Rev.D 77, 083515 

(2008). 

[8] R. J. Scherrer and A. A. Sen, Phys.Rev.D 78, 067303 
(2008). 

[9] A. Ali, M. Sami and A. A. Sen. larXiv:0904.1070l [astro- 
ph.CO]. 

[10] J. A. Frieman, et al., Phys.Rev.Lett. 7 5, 2077 (1995). 
[11] A. Abrahamse, et al.. larXiv:0712.2879l 
[12] U. Alam, V. Sahni and A. A. Staronbinsky, JCAP 0304, 
002 (2003). 

[13] A. Sen, JHEP 9910, 008 (2002); A. Sen, 
Mod.Phys.Lett.A 17, 1797 (2002); M. R. Garousi, 
Nucl. Phys. B584, 284 (2000); E. A. Bergshoeff, M. de 
Roo, T. C. de Wit, E. Eyras, S. Panda, JHEP 0005, 



009 (2000); J. Kluson, Phys.Rev. D 62, 126003 (2000); 
D. Kutasov and V. Niarchos, Nucl.Phys. B 666, 56 
(2003). 

[14] E. J. Copeland, M. R. Garousi, M. Sami , S. Tsu- 
jikawa, Phys. Rev. D71, 043003 (2005); S. Tsu- 
jikawa and M. Sami, Phys. Lett. B603, 113(2004); M. 
Sami, N. Savchenko and A. Toporensky, Phys. Rev. 
D70, 123528(2004); G. N. Felder, Lev Kofman and 
A. Starobinsky,JHEP 0209,026(2002); J. S. Bagla, 
H. K. Jassal and T. Padmanabhan, Phys. Rev. D 67, 
063504 (2003); L. R. W. Abramo and F. Finelli, Phys. 
Lett. B 575, 165 (2003); J. M. Aguirregabiria and 
R. Lazkoz, Phys. Rev. D 69, 123502 (2004); 

[15] J. Simon, L. Verde, R. Jimenez, Phys.Rev.D 71, 123001 
(2005). 

[16] L. Samushia and B. Ratra, Ap.J. 650, L5 (2006). 
[17] M. Hicken et al., Astrophys. J. 700 1097 (2009). 
[18] D. J. Eisensetin, et al., Ap.J. 633, 560 (2005). 
[19] W. J. Percival, et al., MNRAS 381, 1053 (2007). 
[20] J. Kratochvil, A. Linde, E. Linder and M. Shmakova, 
JCAP 0407, 001 (2004). 



